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SIR epidemics in monogamous populations with recombination 

Damian H. Zanett^ 

We study the propagation of an SIR (susceptible-infectious-recovered) disease over an 
agent population which, at any instant, is fully divided into couples of agents. Couples 
are occasionally allowed to exchange their members. This process of couple recombination 
can compensate the instantaneous disconnection of the interaction pattern and thus allow 
for the propagation of the infection. We study the incidence of the disease as a function of 
its infectivity and of the recombination rate of couples, thus characterizing the interplay 
between the epidemic dynamics and the evolution of the population's interaction pattern. 



I. Introduction 

Models of disease propagation are widely used to 
provide a stylized picture of the basic mechanisms 
at work during epidemic outbreaks and infection 
spreading 1.. Within interdisciplinary physics, 
they have the additional interest of being closely re- 
lated to the mathematical representation of such di- 
verse phenomena as fire propagation, signal trans- 
mission in neuronal axons, and oscillatory chemical 
reactions [2] . Because this kind of model describes 
the joint dynamics of large populations of interact- 
ing active elements or agents, its most interesting 
outcome is the emergence of self-organization. The 
appearance of endemic states, with a stable finite 
portion of the population actively transmitting an 
infection, is a typical form of self-organization in 
epidemiological models [5]. 

Occurrence of self-organized collective behavior 
has, however, the sine qua non condition that in- 
formation about the individual state of agents must 
be exchanged between each other. In turn, this re- 
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quires the interaction pattern between agents not 
to be disconnected. Fulfilment of such requirement 
is usually assumed to be granted. However, it is 
not difficult to think of simple scenarios where it 
is not guaranteed. In the specific context of epi- 
demics, for instance, a sexually transmitted infec- 
tion never propagates in a population where sex- 
ual partnership is confined within stable couples or 
small groups [I]. 

In this paper, we consider an SIR (susceptible- 
infectious-recovered) epidemiological model [3] in 
a monogamous population where, at any instant, 
each agent has exactly one partner or neighbor 
[HE]- The population is thus divided into cou- 
ples, and is therefore highly disconnected. How- 
ever, couples can occasionally break up and their 
members can then be exchanged with those of other 
broken couples. As was recently demonstrated for 
SIS models [S1E], this process of couple recombi- 
nation can compensate to a certain extent the in- 
stantaneous lack of connectivity of the population's 
interaction pattern, and possibly allow for the prop- 
agation of the otherwise confined disease. Our main 
aim here is to characterize this interplay between 
recombination and propagation for SIR epidemics. 

In the next section, we review the SIR model and 
its mean field dynamics. Analytical results are then 
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provided for recombining monogamous populations 
in the limits of zero and infinitely large recombina- 
tion rate, while the case of intermediate rates is 
studied numerically. Attention is focused on the 
disease incidence -namely, the portion of the pop- 
ulation that has been infectious sometime during 
the epidemic process- and its dependence on the 
disease infectivity and the recombination rates, as 
well as on the initial number of infectious agents. 
Our results are inscribed in the broader context of 
epidemics propagation on populations with evolv- 
ing interaction patterns [HISllMTTj. 



II. SIR dynamics and mean field de- 
scription 

In the SIR model, a disease propagates over a pop- 
ulation each of whose members can be, at any given 
time, in one of three epidemiological states: suscep- 
tible (S), infectious (I), or recovered (R). Suscep- 
tible agents become infectious by contagion from 
infectious neighbors, with probability A per neigh- 
bor per time unit. Infectious agents, in turn, be- 
come recovered spontaneously, with probability 7 
per time unit. The disease process S — > I — > R 
ends there, since recovered agents cannot be in- 
fected again [3]- 

With a given initial fraction of S and I-agents, 
the disease first propagates by contagion but later 
declines due to recovery. The population ends in 
an absorbing state where the infection has disap- 
peared, and each agent is either recovered or still 
susceptible. In this respect, SIR epidemics dif- 
fers from the SIS and SIRS models, where -due 
to the cyclic nature of the disease,- the infection 
can asymptotically reach an endemic state, with a 
constant fraction of infectious agents permanently 
present in the population. 

Another distinctive equilibrium property of SIR 
epidemics is that the final state depends on the 
initial condition. In other words, the SIR model 
possesses infinitely many equilibria parameterized 
by the initial states. 

In a mean field description, it is assumed that 
each agent is exposed to the average epidemiolog- 
ical state of the whole population. Calling x and 
y the respective fractions of S and I-agents, the 
mean field evolution of the disease is governed by 



the equations 

x = -k\xy, , 
y = k\xy - y, ' ' 

where k is the average number of neighbors per 
agent. Since the population is assumed to re- 
main constant in size, the fraction of R-agents is 
2=1 — x — y. In the second equation of Eqs. 
(fTJ), we have assigned the recovery frequency the 
value 7 = 1, thus fixing the time unit equal to 7 -1 , 
the average duration of the infectious state. The 
contagion frequency A is accordingly normalized: 
A/7 — > A. This choice for 7 will be maintained 
throughout the remaining of the paper. 

The solution to Eqs. Q implies that, from an 
initial condition without R-agents, the final frac- 
tion of S-agents, x* , is related to the initial fraction 
of I-agents, yo, as [I] 

x* = l-(k\)- 1 \og[(l-y )/x*}. (2) 

Note that the final fraction of R-agents, z* = 1 — 
x* , gives the total fraction of agents who have been 
infectious sometime during the epidemic process. 
Thus, z* directly measures the total incidence of 
the disease. 

The incidence z* as a function of the infectiv- 
ity fcA, obtained from Eq. ^ through the stan- 
dard Newton-Raphson method for several values 
yo of the initial fraction of I-agents, is shown in 
the upper panel of Fig. [TJ As expected, the dis- 
ease incidence grows both with the infectivity and 
with j/o- Note that, on the one hand, this growth 
is smooth for finite positive yo. On the other hand, 
for y (but y ^ 0) there is a transcritical 
bifurcation at kX = 1. For lower infectivities, the 
disease is not able to propagate and, consequently, 
its incidence is identically equal to zero. For larger 
infectivities, even when the initial fraction of I— 
agents is vanishingly small, the disease propagates 
and the incidence turns out to be positive. Finally, 
for yo — no agents are initially infectious, no in- 
fection spreads, and the incidence thus vanishes all 
over parameter space. 

III. Monogamous populations with 
couple recombination 

Suppose now that, at any given time, each agent 
in the population has exactly just one neighbor or, 
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Figure 1: SIR epidemics incidence (measured by 
the final fraction of recovered agents z*) as a func- 
tion of the infectivity (measured by the product 
of the mean number of neighbors times the infec- 
tion probability per time unit per infected neigh- 
bor, kX), for different initial fractions of infec- 
tious agents, yo. Upper panel: For the mean field 
equations (UJ. Lower panel: For a static (non- 
recombining) monogamous population, described 
by Eqs. © with r = 0. 



in other words, that the whole population is al- 
ways divided into couples. In reference to sexually 
transmitted diseases, this pattern of contacts be- 
tween agents defines a monogamous population [5] . 
If each couple is everlasting, so that neighbors do 
not change with time, the disease incidence should 
be heavily limited by the impossibility of propagat- 
ing too far from the initially infectious agents. At 
most, some of the initially susceptible agents with 
infectious neighbors will become themselves infec- 
tious, but spontaneous recovery will soon prevail 
and the disease will disappear. 



If, on the other hand, the population remains 
monogamous but neighbors are occasionally al- 
lowed to change, any I-agent may transmit the 
disease several times before recovering. If such 
changes are frequent enough, the disease could per- 
haps reach an incidence similar to that predicted by 
the mean field description, Eq. JT]) (for k = 1, i.e. 
with an average of one neighbor per agent). 

We model neighbor changes by a process of cou- 
ple recombination where, at each event, two cou- 
ples and (m, n) are chosen at random and 
their partners are exchanged [HIE]- The two pos- 
sible outcomes of recombination, either (i, m) and 
(j, n) or (i,n) and (J,m), occur with equal prob- 
ability. To quantify recombination, we define r as 
the probability per unit time that any given couple 
becomes involved in such an event. 

A suitable description of SIR epidemics in 
monogamous populations with recombination is 
achieved in terms of the fractions of couples of 
different kinds, mss, msi, mil, mm,, wrr, and 
msR = 1 — msi — mu — m m — m,RR- Evolution 
equations for these fractions are obtained by con- 
sidering the possible transitions between kinds of 
couples due to recombination and epidemic events 
[7]. For instance, partner exchange between two 
couples (S,S) and (I,R) which gives rise to (S,I) 
and (S,R), contributes positive terms to the time 
derivative of msi and msR, and negative terms to 
those of mss and mm, all of them proportional to 
the product mssmiR- Meanwhile, for example, con- 
tagion can transform an (S,I)-couple into an (1,1)- 
couple, with negative and positive contributions to 
the variations of the respective fractions, both pro- 
portional to msi- 

The equations resulting from these arguments 
read 

m S s = rA sm , 

msi = rB sm - (1 + A)m S i, 

mu = rA ms + \m S i - 2m u , . . 

mm = rBiRS + 2m n - mm, ° 

m RR = rA RS i + mm, 

msR = rB RS i + msi- 

For brevity, we have here denoted the contribution 
of recombination by means of the symbols 

A+jh = (mi J +m ih ) 2 /4-m ii (m :j: j+m : j h +m hh ), (4) 
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and 

Bijh = (2mu + m lh )(2m JJ + m. ]h )/2 

-rriij (rriij + m ih + m jh + m hh )/2, (5) 

with i, j, h e {S, I, R}. The remaining terms stand 
for the epidemic events. In terms of the couple 
fractions, the fractions of S, I and R-agents are 
expressed as 

x = m ss + (msi + m SR )/2, 

y = mn + (msi + m m )/2, (6) 

z = m RR + (m SR + m m )/2. 

Assuming that the agents with different epidemi- 
ological states are initially distributed at random 
over the pattern of couples, the initial fraction of 
each kind of couple is mss(0) = x§, msi(0) = 
2x y , mn(0) = y%, m m (0) = 2y z , m RR (0) = z%, 
and msR,(0) = 2xqZq, where xq, yo and Zq are the 
initial fractions of each kind of agent. 

It is important to realize that the mean field-like 
Eqs. (J3I) to ((HI) are exact for infinitely large popula- 
tions. In fact, first, pairs of couples are selected at 
random for recombination. Second, any epidemic 
event that changes the state of an agent modifies 
the kind of the corresponding couple, but does not 
affect any other couple. Therefore, no correlations 
are created by either process. 

In the limit without recombination, r = 0, the 
pattern of couples is static. Equations §5$ become 
linear and can be analytically solved. For asymp- 
totically long times, the solution provides -from the 
third of Eqs. ([5])- the disease incidence as a func- 
tion of the initial condition. If no R-agents are 
present in the initial state, the incidence is 

z* = (l + A)- 1 [l + A(2-y )]2/o. (7) 

This is plotted in the lower panel of Fig. Q] as a 
function of the infectivity kX = A, for various val- 
ues of the initial fraction of I-agents, yo. When 
recombination is suppressed, as expected, the in- 
cidence is limited even for large infectivities, since 
disease propagation can only occur to susceptible 
agents initially connected to infectious neighbors. 
Comparison with the upper panel makes apparent 
substantial quantitative differences with the mean 
field description, especially for small initial frac- 
tions of I-agents. 



Another situation that can be treated analyti- 
cally is the limit of infinitely frequent recombina- 
tion, r — > oo. In this limit, over a sufficiently short 
time interval, the epidemiological state of all agents 
is virtually "frozen" while the pattern of couples 
tests all possible combinations of agent pairs. Con- 
sequently, at each moment, the fraction of couples 
of each kind is completely determined by the in- 
stantaneous fraction of each kind of agent, namely, 

m S s = x 2 , m S i = 2xy, m n = y 2 , , 
mm = 2yz, m RR = z 2 , m SR = 2xz. 

These relations are, of course, the same as quoted 
above for uncorrelated initial conditions. 

Replacing Eqs. ([8]) into ([3]) we verify, first, that 
the operators Aijh and B^h vanish identically. The 
remaining of the equations, corresponding to the 
contribution of epidemic events, become equivalent 
to the mean field equations ([1]). Therefore, if the 
distributions of couples and epidemiological states 
are initially uncorrelated, the evolution of the frac- 
tion of couples of each kind is exactly determined 
by the mean field description for the fraction each 
kind of agent, through the relations given in Eqs. 
©• 

For intermediate values of the recombination 
rate, < r < oo, we expect to obtain incidence lev- 
els that interpolate between the results presented 
in the two panels of Fig. [TJ However, these can- 
not be obtained analytically. We thus resort to the 
numerical solution of Eqs. J3|). 

IV. Numerical results for recombin- 
ing couples 

We solve Eqs. ([3]) by means of a standard fourth- 
order Runge-Kutta algorithm. The initial condi- 
tions are as in the preceding section, representing 
no R-agents and a fraction j/o of I-agents. The dis- 
ease incidence z* is estimated from the third equa- 
tion of Eqs. ((6|), using the long-time numerical so- 
lutions for m RR , msR, and mi R . In the range of 
parameters considered here, numerical integration 
up to time t — 1000 was enough to get a satisfac- 
tory approach to asymptotic values. 

Figure [5] shows the incidence as a function of in- 
fectivity for three values of the initial fraction of 
I-agents, y — > 0, yo — 0.2 and 0.6, and several 
values of the recombination rate r. Numerically, 
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Figure 2: SIR epidemics incidence as a function 
of the infectivity for three initial fractions of infec- 
tious agents, yo, and several recombination rates, 
r. Mean field (m. f.) results are also shown. The 
insert in the upper panel displays the boundary be- 
tween the phases of no incidence and positive inci- 
dence for yo — » 0, in the parameter plane of infec- 
tivity vs. recombination rate. 



positive values of r give rise to incidences between 
those obtained for a static couple pattern (r = 0) 
and for the mean field description. Note that sub- 
stantial departure from the limit of static couples 
is only got for relatively large recombination rates, 
r > 1, when at least one recombination per cou- 
ple occurs in the typical time of recovery from the 
infection. 

Among these results, the most interesting situa- 
tion is that of a vanishingly small initial fraction of 
I-agents, yo — > 0. Figure [3] shows, in this case, the 
epidemics incidence as a function of the recombi- 
nation rate for several fixed infectivities. We recall 
that, for yo —5" 0, the mean field description predicts 
a transcritical bifurcation between zero and posi- 
tive incidence at a critical infectivity X = 1, while 
in the absence of recombination the incidence is 
identically zero for all infectivities. Our numerical 
calculations show that, for sufficiently large values 
of r, the transition is still present, but the criti- 
cal point depends on the recombination rate. As 
r grows to infinity, the critical infectivity decreases 
approaching unity. 
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Figure 3: SIR epidemics incidence as a function 
of the recombination rate r for a vanishingly small 
fraction of infectious agents, yo - > 0, and several 
infectivities A. 



the limit yo — > has been represented by taking 
yo = 10 -9 . Within the plot resolution, smaller val- 
ues of yo give identical results. Mean field (m. f.) 
results are also shown. As expected from the an- 
alytical results presented in the preceding section, 



Straightforward linearization analysis of Eqs. © 
shows that the state of zero incidence becomes un- 
stable above the critical infectivity 
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This value is in excellent agreement with the nu- 
merical determination of the transition point. Note 
also that Eq. © predicts a divergent critical infcc- 
tivity for a recombination rate r = 1. This implies 
that, for < r < 1, the transition is absent and 
the disease has no incidence irrespectively of the 
infectivity level. For yo — > 0, thus, the recombina- 
tion rate must overcome the critical value r c = 1 
to find positive incidence for sufficiently large infec- 
tivity. The critical line between zero and positive 
incidence in the parameter plane of infectivity vs. 
recombination rate, given by Eq. (j^]), is plotted in 
the insert of the upper panel of Fig. [5J 

V. Conclusions 

We have studied the dynamics of SIR epidemics in 
a population where, at any time, each agent forms a 
couple with exactly one neighbor, but neighbors are 
randomly exchanged at a fixed rate. As it had al- 
ready been shown for the SIS epidemiological model 
[HE] , this recombination of couples can, to some de- 
gree, compensate the high disconnection of the in- 
stantaneous interaction pattern, and thus allow for 
the propagation of the disease over a finite portion 
of the population. The interest of a separate study 
of SIR epidemics is based on its peculiar dynamical 
features: in contrast with SIS epidemics, it admits 
infinitely many absorbing equilibrium states. As 
a consequence, the disease incidence depends not 
only on the infectivity and the recombination rate, 
but also on the initial fraction of infectious agents 
in the population. 

Due to the random nature of recombination, 
mean field-like arguments provide exact equations 
for the evolution of couples formed by agents in 
every possible epidemiological state. These equa- 
tions can be analytically studied in the limits of 
zero and infinitely large recombination rates. The 
latter case, in particular, coincides with the stan- 
dard mean field description of SIR epidemics. 

Numerical solutions for intermediate recombina- 
tion rates smoothly interpolate between the two 
limits, except when the initial fraction of infectious 
agents is vanishingly small. For this special situ- 
ation, if the recombination rate is below one re- 
combination event per couple per time unit (which 
equals the mean recovery time), the disease does 
not propagate and its incidence is thus equal to 



zero. Above that critical value, a transition ap- 
pears as the disease infectivity changes: for small 
infectivities the incidence is still zero, while it be- 
comes positive for large infectivities. The critical 
transition point shifts to lower infectivities as the 
recombination rate grows. 

It is worth mentioning that a similar transition 
between a state with no disease and an endemic 
state with a permanent infection level occurs in SIS 
epidemics with a vanishingly small fraction of in- 
fectious agents [H[2]- For this latter model, how- 
ever, the transition is present for any positive re- 
combination rate. For SIR epidemics, on the other 
hand, the recombination rate must overcome a crit- 
ical value for the disease to spread, even at very 
large infectivities. 

While both the (monogamous) structure and the 
(recombination) dynamics of the interaction pat- 
tern considered here are too artificial to play a role 
in the description of real systems, they correspond 
to significant limits of more realistic situations. 
First, the monogamous population represents the 
highest possible lack of connectivity in the interac- 
tion pattern (if isolated agents are excluded). Sec- 
ond, random couple recombination preserves the 
instantaneous structure of interactions and does 
not introduce correlations between the individual 
epidemiological state of agents. As was already 
demonstrated for SIS epidemics and chaotic syn- 
chronization [7], they have the additional advan- 
tage of being analytically tractable to a large ex- 
tent. Therefore, this kind of assumption promises 
to become a useful tool in the study of dynamical 
processes on evolving networks. 
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